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Exact analytical results (in terms of Bessel functions) for the bandgaps, reflectance, and trans- 
mittance of one-dimensional photonic crystals with a sawtooth refractive index profile on the period 
are derived for the first time. This extends a group of exactly solvable models of periodic refractive 
indices. The asymptotic approximations of the above exact results have been also obtained. 
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I. INTRODUCTION 

Recently there has been a renewal of interest in the 
properties of ID photonic crystals based on graded in- 
dex slabs. Rauh et al. [J have considered in some detail 
optical properties of periodic systems consisting of slabs 
whose permittivities n 2 (z) increase linearly with depth z. 
They also considered periodic systems in which alternat- 
ing layers have quadratically increasing and decreasing 
permittivities. Their work pointed out the relevance of 
such ID photonic crystals for technological applications. 
They stressed in particular that to a good approximation 
all the bandgaps have the same width. 

In this paper we consider the properties of ID photonic 
crystal constructed of layers with linearly increasing re- 
fractive index n(z), giving a quadratic increase of the 
permittivity n 2 (z). To be precise, for a slab extending 
from z = to z = d we take 

n(z) = n a + (n b - n a )- = nb — - (z - z ), (1) 
a a 

with zq = —n a d/ (nt,—n a ). As seen in Fig. 1 the refractive 
index is periodic, increasing linearly from n a to rib inside 
each layer, and falls sharply at the cell boundaries. From 



n(z), 



/ib 





d 2d 



(N-1)d Nd 



data in Adachi Q one finds that the permittivity of the 
ternary alloy Al^Gai-^As is approximately linear in the 
relative content x > 0.4 of Al for wavelengths 10000 A 
> A > 6000 A, so photonic crystals with linearly graded 
refractive index layers are within present technological 
possibilities. 

We will express the dispersion equation for the 
bandgaps of these photonic crystals as well as their re- 
flectance and transmittance in terms of Bessel functions 
of fractional order. This adds to a group of exactly 
solvable models of periodic refractive indices which is 
restricted at the moment to periodic structures with a 
step-layered profile of the index on the period (binary 
and ternary photonic crystals) and by some specific two- 
parameter sinusoidal periodic potentials 0, H[ . By means 
of asymptotic expansions we derive approximate sinu- 
soidal expressions showing that the bandgaps of our saw- 
tooth periodic potential tend to a constant width as the 
wavelength decreases. For simplicity, as in Rauh's work 
[l| , we restrict our attention to normal propagation. 



II. EXACT RESULTS 

Description of optical wave propagation of linearly po- 
larized light of circular frequency oj through a periodic 
structure, with real refractive index n(z + d) = n(z), in 
case of normal incidence, reduces from Maxwell's equa- 
tions to the Hclmholtz equation: 



e>> 



k 2 n 2 (z)E(z) = 0, 



(2) 



where the total electric and magnetic fields of propagat- 
ing light arc expressed as 



FIG. 1: ID photonic crystal with a sawtooth refractive index. 
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(3) 
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Applying the transfer matrix method for periodic poten- 
tials 0, the field E(z) and its derivative dE(z)/dz at the 
edge points of the system z = and z = Nd are related 
by 



E {Nd) 
E' (Nd) 



v' (d) -v (d) 
-u' (d) u(d) 



N 



25(0) 
E'(0) 



(4) 



where u(z) and v(z) are the so-called normalized solu- 
tions, which satisfy the boundary conditions 



(5) 



The TV'th power W N of a unimodular 1^-matrix as oc- 
curs in Eq. (4) can be expressed in terms of W and the 
unit matrix [1] as shown in Ref. 
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[1], 



(6) 



where the dispersion equation for the Bloch phase <p takes 
the form 



cos <f) = i Tr W = - [u(d) + v'(d)] 



(7) 



The advantage of the transfer matrix method is that once 
we have solved for W of a single layer, the result for N 
layers is trivial, by use of Eq. (6). 

In dimensionless units, Eq. (2) with the refractive in- 
dex defined by Eq. (1), takes the form 



d 2 E(S) 1 2 / \ 

where the dimensionless variable z is 

,-, v 2k (n b -n a ) 
z = yH^- z ), 7= 



(8) 



(9) 



It is straightforward to prove (see Appendix A) that the 
normalized solutions of Eq. (8) are Bessel functions of 
order v = ±1/4, multiplied by z 1 ! 2 , i.e. 



u{z) 



v{z) 



2 1 /4r(l/4) 



2 3 /4r(3/4) 



/4 



/i 



48 

£! 

80 " 



(10) 



Using standard relations among Bessel functions, see 
chapter 9 in Q, one finds 



du(z) 

dz, 
dv(z) 

dz 



1 g3/2 j Oi 

2 5 / 4 r(l/4) 3/4 V 4 

1 =3/2 j 

2 7 /4r(3/4)" 3/4 l 4 



(11) 



If we now change the variable z in Eq. (10) back to the 
original variable z, we obtain another set of the funda- 
mental solutions, u(z) and v(z), of Eq. (2) 



u(z) 



v(z) 



2!/4r(i/4) 



2 3 /4r(3/4) 



Wl{z - Zo)} 1/2 J-l/4 



Wl{* - Zo)} 1/2 Jl/4 



l(z - z a f 



(12) 



To construct the matrix W of Eq. (4), one has to ex- 
press u(z) and v(z) in terms of u(z) and v(z). Taking 
into account Eq. (5), one obtains 



5'(0) . fi'(0) , 
u \ z ) = ~/n\ u ( z > z 77v\ v \ z )' 



w(0) 



w(0) 



v(0) „, s 6(0) , 



(13) 



w(0) 



w(0) 



where the Wronskian w(z) of the fundamental solutions 
u(z) and v(z) of Eq. (2) at the point z = is 



w(0) = u{0)v'(0) - v(0)u'(0) = V7- (14) 

With the aid of newly introduced dimensionless parame- 
ters, 



H z 



2k Q dn 2 



Zb = Vl(d- z ) 



n b - n a 
odnl 



2kodr 2 ^ 1/2 



n b - n a 



(15) 



the values of u(0), v(0), u'(0), v'(0) take the form 

2(0) 
0(0) 



" ,1/2 j 

2V4r(l/4) a 1/4 V 4 



" 1/2 T [ 

23/4r(3/4) 11 1/4 I 4 



fi'(0) 



_ du(z) 



v'(0) 



dz 
dv(z) 



= -V7 



n 3/2 t / z " 



dz 



z=0 



z=Q 



2 5 / 4 r(i/4) 



z a J 3/4 



i— " 3/2 t I 

^ 7 2V4r( 3 /4) Za J - 3/4 lvT 



(16) 



while at z = d. u(d), v(d), u'(d), v'(d) are given by 



u{d) 
v(d) 



* 1/2 j 

2 i/4r(i/ 4 ) Zh J - 1/4 



K 1/2 . 

23/4r(3/4) Zfc 1/4 V 4 



M '(d) EE *1 



-V7 
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^ 2V4r( 3 /4) Zb3/2j - 3/4 



(17) 
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The elements of the W^-matrix arc then 



W u = v'{d) = 
Wai 



2 5 /2 



Z a Zb 



J 1/4 (4/4) J 3 /4 (4/4) + J-1/4 (4/4) J-3/4 (4/4) 



u'(d) 



fc 

k Wi2 

w 22 



fc 2 5 / 2 



^75- y/n a n b z a z b 



2 a ' z y rib 



J_ 3 /4 (4/4) J3/4 (4/4) - J3/4 (4/4) J-3/4 (4/4) 

ZaZb Jl/4 (4/4) J-1/4 (4/4) - J-1/4 (4/4) Jl/4 (4/4) 

J-3/4 (4/4) J-1/4 (4/4) + J3/4 (4/4) J1/4 (4/4) 



(18) 



while the dispersion equation (7) takes the form 



2 7 /2 



J-i J-3/4 (4/4) J-1/4 (4/4) + J— J3/4 (4/4) Jl/4 (4/4) 
71 5 y 715 

+ Ji/4 (4/4) J3/4 (4/4) + M J-1/4 (4/4) J-3/4 (4/4) 



(19) 



Eq. [T5] is the exact dispersion relation for our sawtooth ID photonic crystal. The amplitude reflection r and trans- 
mission t coefficients arc 



W N_'J^ W N +i 



'ex TTfA^ 
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(20) 



W N+'J^ W N + i 
n in 



k n cx Wr 2 - 



k Tlh 



-W, 



N 



where the elements of the W-matrix are given in exact form by Eq. (18) and the elements of the VF^-matrix are 
obtained by use of Eq. (6). The indices of refraction in the incident and exit media are n ln and n ex as seen in Fig. 1. 



III. ASYMPTOTIC ESTIMATES 



The motivation for exploring those is to have simple analytic approximations for the Bloch phase cos <f>, as well as 
for the reflection and transmission coefficients. Using Hankel's asymptotic formula, section 9.2 of Ref. Q, one obtains 
the elements of the W matrix in the form 



, rib 
n c 

n. 
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ko 



k Wi2 = 
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1 

\Jn a n b 
fn_ 
m, 
fn 
rib 



P 1/4 (4/4) P 3 /4(4/4) + Ql/4(4/4) 03/4(4/4) 

Pi/4(4/4) Q 3/4 (4/4) - 0i/ 4 (4/4) P 3 /4(4/4) 

P 3/4 (4/4) ^3/4(4/4) + 03/4(4/4) 03/4(4/4) 
P 3/4 (4/4) 03/4(4/4) - 03/4(4/4) ^3/4(4/4) 

P 1/4 (4/4) A/ 4 (4/4) + Qi/ 4 (4/4) 0i/ 4 (4/4) 

a /4 (4/4) 0i /4 (4/4) - a /4 (4/4) 0i /4 (4/4) 

P 1/4 (z fc 2 /4) P 3/4 (4/4) + 0i/ 4 (4/4) 03/4(4/4) 
P 3 /4(4/4) 0i/ 4 (4/4) - Pi/ 4 (4/4) 03/4(4/4) 



cos(z 2 /4-z 2 /4) 
sin (z 2 /4 - zl/A) , 
sin (z 2 /4 - z 2 /4) 
cos (z 2 /4 - z 2 /4) , 
sin(z 2 /4-z 2 /4) 
cos (4/4 - z b 2 /4) , 
cos (z 2 /4-z b 2 /4) 



sin 



(4/4-4/4) 



(21) 
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where 



Pu{z) = 1 

Q v {z) 



{Ay 2 - \){Ay 2 - 9) {Ay 2 - l){Av 2 - 9){Ay 2 - 25)(4^ 2 - 49) 
2!(8z) 2 4!(8z) 4 
4^ 2 -l (4^-1)(4i/ 2 -9)(4j/ 2 -25) 



8.2 



3!(8z)= 



and the dispersion equation in the form 



cos <f) = - > 


I n a 


P 1/4 (z 2 /4)P 3/4 (z 2 /4) 


+ Qi/ 4 (4/4) Q 3/4 (4/4) 


cos 


(4/4 " 


- 4/4) 


1 

+ 2\ 


' n b 


P 1/4 (z 2 /4)P 3/4 (z 2 /4)- 


^Qi/4(z 2 /4)Q 3/4 (z 2 a /4) 


cos 


(4/4 - 


4/4) 




[rib 
' n a 


^i/ 4 (4/4) Q 3 /a(4/*) 


-Q 1/4 (z 2 /4)P 3/4 (z 2 /4) 


sin 


:^/4 - 


4/4) 




fn^ 
' n b 


P 3/4 (4/4)Q 1/4 (z b 2 /4) 


-P 1/4 (z 2 /4)Q 3/4 (z 2 /4) 


sin 


:^/4- 


4/4) 



(22) 



(23) 



Eqs. (21) - (23) are the expressions we sought. By truncating the terms P„ and Q v various asymptotic approxima- 
tions are obtained. In particular, after some trigonometric simplifications, the matrix elements become in the first 
approximation (P = 1, Q = 0) 



Wax 
k 

k W 12 
W 22 



cos (z 2 /4 - z 2 /4) , 

ria 

-V n a n-b sin (4/4 - 4/ 4 ) 
sin (z 2 /4 - z 2 /4) , 



fTi 



^cos(z 2 /4-z 2 /4), 

Tib 



in the second (P = 1, Q = {v 2 - 0.25)/2z ), 



Wu = 

W21 = 
k 

k W 12 = 



(z 2 a /A - z 2 /A) + (_ + _) sin (z 2 /4 - z 2 /4) 



IT" 



while in the third, 



n b 



Wu = J — 



1 



22 



105 



-\Jn a n b 
1 

\Jri a rib 

n_ 

n>b 



in (z 2 /4 - z 2 /4) +(_-_) cos (z 2 /4 - z 2 /4) 



cos 



in (z 2 /4 - z 2 /4) -^—-—J cos (z 2 /4 - z 2 /4) 

(4/4 - 4/4) - + 4) sin {zl/A " z&2/4) 



135 



15 



128z 4 128z 4 64z 2 z 2 



(z 2 /4 - 2fc 2 /4) + (_ + _) sin (z 2 /4 - z 2 /4) 



W 2 



k 
k Wi 2 



II 



22 



1 

Fn 
n h 



135 



135 



25 



128z 4 128z 4 64z 2 z 2 



105 



105 



128z 4 128z 4 64z 2 z 2 



in (z 2 /4 - z 2 /4) + (___) cos (z 2 /4 - z 2 /4) 
sin (z 2 /4 - z 2 /4) - ( - ^ ) cos (z 2 /4 - z 2 /4) 



1 



135 



105 



15 



z 2 



COS 



;4/4 - 4/4) 



J. 



\8z 2 8z 2 



sin 



;4/4-4/4) 



(24) 



(25) 



(26) 



As one can see, the dispersion equation cos <j> = i (W11+W22) takes a particularly simple form in the first approximation 

cos <fi = ^ (7? + cos (z 2 /4 - z 2 /4) . (27) 



5 



Figs. 2, 3, 4 illustrate the above approximations. In 
the second and third approximations, cos0 of Eq. (23) 
takes the form 



cos 



02,3 = A 2 . 3 cos(z*/4 - zl/A + 6 2 , 3 ), (28) 



which allows for a shift of both the band centers and 
their widths. For our choice of n&/n a = 2, the amplitude 
in the first approximation A\ = 1.06 accounts for the 
narrow band gaps observed in the drawings. 




FIG. 2: (colour online) Exact results (grey lines), see Eqs. (18- 
20), vs. the first asymptotic approximation (red lines), see 
Eq. (24), for the Bloch phase, reflectance, and transmittance, 
for a sawtooth periodic structure with n a = 1.5, nt, = 3, and 
the number of periods N — 6, in the case of normal incidence. 
The refractive indices outside the structure are ru n = 1.0 for 
2 < 0, and n ox = 1.0 for z > Nd. 



In Appendix B we also compare the exact results to 
approximations provided by binary photonic crystals (re- 
placing the ramp by discrete steps). 




FIG. 3: (colour online) Exact results (grey lines), see Eqs. (18- 
20), vs. the second asymptotic approximation (green lines), 
see Eq. (25), for the Bloch phase, reflectance, and transmit- 
tance, for a sawtooth periodic structure with parameters as 
in Fig. 2. 



IV. BANDGAP ANALYSIS 

Fig. [U shows that the dispersion equation for cos </> in 
the first approximation is already in quite good agree- 
ment with the exact results. We now discuss some 
of the predictions for the bandgaps that follow from 
it. The bandgap edges are defined by the condition 
cos<f> = ( — l) f/ , where the index q = 1,2,... numbers 
the bandgaps. This condition applied to Eq. (27) leads 
to expressions for the right and left fcg bandedges in 
the first approximation: 



k r d 



qir 



k l d 



qir 



(29) 
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FIG. 4: (colour online) Exact results (grey lines), see Eqs. (18- 
20), vs. the third asymptotic approximation (blue lines), see 
Eq. (26), for the Bloch phase, reflectance, and transmittance, 
for a sawtooth periodic structure with parameters as in Fig. 2. 



is inherent in use of an asymptotic series), all three ap- 
proximations provide one with quite accurate estimates 
and the second one balances accuracy and simplicity. 



0.008- 
0.006- 
0.004 - 
0.002- 
0.000- 
-0.002- 
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FIG. 5: (colour online) Dimensionless difference, (fe§,a ~~ &o) d, 
a = 1, 2, 3, between exact results, kg d, and a given asymptotic 
approximation, fcg,a d, for bandgap centers of a sawtooth peri- 
odic array with refractive index parameters n a = 1.5, n;, = 3.0 
vs. the bandgap index q for the first eight bandgaps. From 
top to bottom the lines are — ^o) d ( re d), (&§ — kg)d 
(grey), (fco, 3 - ko)d (blue), (fcg j2 - feg)d (green). 



T3 
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where 

( 2 \ n a + n b 

a = arccos - , n av = . (30) 

\\fn a /n b + n b /n a J * 

Therefore, the centers of the bandgaps /eg and their 
widths w are given in the first approximation by 

7r 2a 
k d = q , wd — . (31) 

This is qualitatively similar to the asymptotic results ob- 
tained by Rauh et al. [IJ. In Fig. 5 we compare the 
accuracy of the three approximations for the band gap 
centers, and in Fig. 6 for their widths. As one can see 
from Eqs. (24) - (26), each approximation incorporates 
new corrections of orders l/z% b and their squares respec- 
tively. With the exception of the first bandgap (but that 



FIG. 6: (colour online) Dimensionless bandgap widths, w d, 
for a sawtooth periodic array with refractive index param- 
eters n a = 1.5, rib = 3.0 vs. the bandgap index q for the 
first eight bandgaps. From top to bottom the lines are the 
second asymptotic approximation (green), the first asymp- 
totic approximation (red), exact result (grey), and the third 
asymptotic approximation (blue). 



V. CONCLUSION 

Wc have derived the exact dispersion equation, reflec- 
tion, and transmission coefficients, in terms of Bcsscl 
functions, for our one-dimensional photonic crystal with 
a sawtooth refractive index profile. Using the Hankcl 
expansion, we worked out the first three asymptotic ap- 
proximations to those coefficients. Even the first approx- 
imation, which has a particularly simple analytic form, 
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provides one with very reasonable estimates. Using the 
results of Appendix A, one can easily work out solutions 
for any monomial behaviour (kon(z)) 2 ~ (z/L) n of the 
permittitivity within a single cell of the periodic system. 
They always involve Bessel functions J± m (z) of order 
m = l/(7i + 2). 
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FIG. 7: (colour online) Exact results (grey lines) for a saw- 
tooth periodic potential with parameters as in Fig. 2 vs. bi- 
nary approximation (black lines) with ni = n a , 112 = rib. 




FIG. 8: (colour online) Like Fig. 7 but with m = (n a +3nt,)/4, 
n 2 = (3n a + n b )/A. 



Appendix A: Solutions of the wave equation for 
power law potentials 

While discussing the WKB matching conditions at a 
classical turning point of order n, Schiff || points out 
that when k 2 (z) = (z/L) n , an exact solution of the sta- 
tionary state wave equation is given by 



= i>"(z) + k 2 (z)^(z), 

where 

1>(z) = A(z) J m (0 , A(z) = VM/Hzj , 



a*) 



k(t)dt, 



(Al) 



(A2) 



The task is to show that J m (0 is a generic Bessel function 
of order m, a result given originally by Langer Q. The 
most used case is a linear turning point with n = 1 => 
m = 1/3, where the solution coincides with the Airy 
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function. The general result makes it easy to work out the 
transfer matrix for a ID photonic crystal whose refractive 
index profile has any monomial "sawtooth shape" in each 
layer. To begin, we note that for our assumed k 2 (z) 



that 



k'(z) 

A\z) 
2A(z)A'(z) 
A 3 (z)A"(z) 



2L /2\(«+2)/2 

~ n + 2 \l) 

n / z n(«-2)/2 

= 2L\l) 
= — = 2m z , 

K 

= 2m = const , 

= - [A{z)A'f = -m 2 



(A3) 



From here we work out ip' and i/j" , and substitute into 
the wave equation. On functions of z, the prime means 
d/dz, but on the Bessel function it means d/d£. We have 



ip' = A'(z)J m (0 + A(z)k(z)JUO, 
V" = A"(z)J m (0 + A(z)k 2 (z)J^(0 

+ \2A'(z)k(z) + A(z)k'(z)]j' m (0 



(A4) 



Substituting into the wave equation and removing the 
factor Ak 2 , we find that the coefficient of J^iO is 



2A'k + Ak' 
Ak 2 



(A5) 



see Eq. (|A8|) for details, while the coefficient of J m (£) is 
1 + [A" /Ak 2 ], where the variable portion is 



A" 
Ak 2 



A 3 (z)A" 
A A k 2 



(A') 2 A 2 
A^k 2 



(A6) 



This verifies that the wave equation has been reduced to 
Bessel's equation for the fractional order m. 



•C(0 + ^4(0 + 



m 



J m (0 = . (A7) 



The regular solution J rn (£) will vanish at the origin when 
m > 0, so the irregular solution J_ m (£) is also required 
if a non-zero initial value is in order. In Eq. (|A5[) we used 



A\z). 
2A'k 



i/k, 
^Ak' 



1 



Ak 2 



k 

A 2 k~ + k 2 ~ £ 



2AA' -- 
2AA' 



(A8) 



Appendix B: Binary Approximations. 

We consider a binary photonic crystal composed of 
layers of constant refractive indices m, n 2 , with corre- 
sponding thicknesses d\, d 2 . The elements of a binary 
VF-matrix are well-known 0] and given by 



W 21 



cos(konidi) cos(fco ^2^2) 
n 2 



■ sin(fconidi) sin(fcon2rf2), 



k 
k W 12 



■ n\ sin(fconidi) cos(fcori2^2) 
n 2 co%(k,Qn\d\) smik^n-id-i) , 
■ sa\(k^n\d\) cos(fco 712^2) 



(Bl) 



"1 



f 



cos(fconidi) sin(fcon2^2), 



22 



112 

cos(konidi) cos(fcon2C?2) 

m 



■ sm(konidi) sin(fcon2(i2)- 



"2 



The simplest binary approximation to the sawtooth 
profile discussed above is given by n\ = n a , n 2 = rib, and 
d\ = d 2 = d/2. Fig. [7] compares the exact results for 
the sawtooth potential with parameters as in Fig. 2 to 
that of the simplest binary approximation. As one can 
see there are some similarities, but the bandgap widths 
differ strongly; for example, the third bandgap of the 
sawtooth potential disappears entirely. A better binary 
approximation, see Fig. [51 is obtained if one follows the 
recipe of Kalotas and Lee taking the layer indices 
to be Tii = (n a + 3nt)/4, n 2 = (3n a + nt,)/4, and widths 
d\ = d 2 = d/2. These values give the correct average 
wave number in each layer. 
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